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The branching fractions for the inclusive Cabibbo-favored K*° and Cabibbo-suppressed K*° de- 
cays of D mesons are measured based on a data sample of 33 pb~ x collected at and around the center- 
of-mass energy of 3.773 GeV with the BES-II detector at the BEPC collider. The branching fractions 
for the decays D +(0) -> i?*°(892)X and D° -> K*°(892)X are determined to be BF(D° -> K*°X) = 
(8.7±4.0±1.2)%, BF(D + -> K*°X) = (23.2±4.5±3.0)% and BF(D° -> K*°X) = (2.8±1.2±0.4)%. 
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An upper limit on the branching fraction at 90% C.L. for the decay D + — > K*°(892)X is set to be 
BF(D+ -> K*°X) < 6.6%. 



PACS numbers: 13.20.Fc,13.25.Ft,13.85.Ni,14.40.Lb 

I. INTRODUCTION 

Although D mesons were found 29 years ago 0, 0, 
the study of charm meson decay properties is still an 
interesting field. Measurement of branching fractions 
for the D meson decay modes containing K*°(K*°) in 
the final states can provide useful information about the 
relative strength of the Cabibbo-favored and Cabibbo- 
suppressed D decays. The total branching fractions for 
the exclusive D decay modes containing K*° are summed 
to be BF(D° -> K*°X) = (8.1 ± 0.8)% and BF(D+ -> 
K*°X) = (23.1 ± 2.0)% {X = any particles) with the 
known branching fractions of neutral and charged D 
mesons |3j . The measurement of branching fractions for 
the inclusive Cabibbo-favored decay D — > K*°X and 
Cabibbo-suppressed decay D — » K*°X can not only serve 
as an independent check on the sum of the branching frac- 
tions for the exclusive D decay modes containing K*° me- 
son in the final states, which indicates the need to search 
for new decay modes, but also provides valuable informa- 
tion in understanding the weak decay mechanism. The 
knowledge of the inclusive D meson decay properties will 
also help one to understand B decays. 

This Letter reports the measurement of branching frac- 
tions for the inclusive Cabibbo-favored decay D — > K*°X 
and Cabibbo-suppressed decay D — > K*°X of neutral 
and charged D mesons using a double tag method de- 
scribed in Section III, based on an analysis of about 33 
pb^ 1 of data collected with the upgraded Beijing Spec- 
trometer (BES-II) in e + e~ annihilation at and around 
= 3.773 GeV. 



II. THE BEIJING SPECTROMETER 

BES is a conventional cylindrical magnetic detector 
P| operated at the Beijing Electron Positron Collider 
(BEPC) 0. BESII is the upgraded version of the 
BES detector 0. A 12-layer vertex chamber (VC) sur- 
rounding the beam pipe provides trigger information. 
A forty-layer main drift chamber (MDC) located out- 
side the VC performs trajectory and ionization energy 
loss (dE/dx) measurement with a solid angle coverage 
of 85% of 47T for charged tracks. Momentum resolution 
of <J p /p = 1.7%y/l + p 2 (p in GeV/c) and dE/dx resolu- 
tion of 8.5% for Bhabha scattering electrons are obtained 
for the data taken at y/s — 3.773 GeV. An array of 48 
scintillation counters surrounds the MDC and measures 
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the time of flight (TOF) of charged tracks with a res- 
olution of about 200 ps for the electrons. Surrounding 
the TOF is a 12-radiation-length, lead-gas barrel shower 
counter (BSC) operated in limited streamer mode, which 
measures the energies of electrons and photons over 80% 
of the total solid angle, and has an energy resolution 
of (Te/E = 0.22/ \[E (E in GeV), spatial resolutions of 
(70 = 7.9 mrad and oz = 2.3 cm for the electrons. Out- 
side of the BSC is a solenoidal magnet which provides a 
0.4 T magnetic field in the central tracking region of the 
detector. Three double-layer muon counters instrument 
the magnet flux return, and serve to identify muons with 
momentum greater than 500 MeV/c. They cover 68% of 
the total solid angle with longitudinal (transverse) spa- 
tial resolution of 5 cm (3 cm) . End-cap time-of-flight and 
shower counters extend coverage to the forward and back- 
ward regions. A Monte Carlo package based on GEANT3 
has been developed for BESII detector simulation and 
comparisons with data show that the simulation is gen- 
erally satisfactory Q. 



III. DATA ANALYSIS 

A. Event Selection 

For each event it is required that at least 2 (but no 
more than 10) charged tracks are well reconstructed in 
the MDC with good helix fits. All tracks, save those 
from Kg decays, must originate from the interaction re- 
gion, which requires that for a charged track, the dis- 
tance of closest approach be less than 2 cm in the xy 
plane, and less than 20 cm in the z direction. For the 
7r + and tt~ from the Kg decay, the secondary vertex po- 
sition is required to be less than 8 cm in the xy plane 
and within ±20 cm in the z direction to the primary 
interaction point. In addition, in order to optimize the 
momentum resolution and charged particle identification, 
a geometry cut | cos6 |< 0.85 (9 is the polar angle of the 
track) is applied. For the charged particle mass assign- 
ment, a combined confidence level calculated using the 
dE/dx and TOF measurements is required to be greater 
than 0.1% for the pion hypothesis. For the kaon identi- 
fication, the confidence level for the kaon hypothesis is 
required to be greater than that for the pion hypothesis. 

The 7T° is reconstructed through the decay 7r° — * 77. 
For the 7 from ir° decay, the energy deposited in the BSC 
is required to be greater than 70 MeV; the electromag- 
netic shower is required to start in the first 5 readout lay- 
ers; and the angle between the 7 and the nearest charged 
track is required to be greater than 22°. 
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B. Singly Tagged D° and D~ 

Around the center-of-mass energy of 3.773 GeV, 
the V'(3770) resonance is produced in electron-positron 
(e + e~) annihilation. The ?/>(3770) lies just above open 
charm pair production threshold and decays predomi- 
nantly into DD pairs. If one D meson is fully recon- 
structed (this is called a singly tagged D event), the other 
D meson must exist on the recoil side. Throughout the 
paper, charge conjugation is implied. 

For the analysis, singly tagged D events are recon- 
structed in three hadronic D° decay modes (K + tt~, 
K + ir~ir~~ir + , and K + ir~ir°) and in nine D~ decay modes 
K°K-, K+K-TT-, K°Tr~' 
t-tt°, 



(K +r K t 



A' H 



TV TT 7T 7T 



and 



7T 7T 



The reconstruction for the singly tagged D events is the 
same as that used in the previous works || 0| ■ 

In order to reduce the background and improve the 
momentum resolution, a center-of-mass energy con- 
straint kinematic fit (lC-fit) is imposed on each mKnir 
(m=0,l,2; n=l,2,3,4) combination. For the singly tagged 
D decay modes with a neutral kaon or a neutral pion 
in the daughter particles, an additional constraint kine- 
matic fit for the Kg — > tt + it~ or ir° — > 77 is also per- 
formed. The kinematic fit probability P(x 2 ) is required 
to be greater than 0.1%. If more than one combination 
satisfies the criteria in an event, only the combination 
with the largest P(x 2 ) is retained. 

Fig. and fig. |2] show the invariant mass spectra for 
mKnir combinations in the singly tagged D° and D~ 
decay modes, which are calculated based on the track 
momentum vectors from the kinematic fit. A maximum 
likelihood fit to the mass spectrum with a Gaussian func- 
tion for the D signal and a special background function 
(ARGUS background shape multiplied by a polynomial 
function) Q to describe backgrounds yields a total num- 
ber of 7033 ± 193 ± 316 singly tagged D° events and 
5321 ± 149 ± 160 singly tagged D~ events (as well as the 
numbers for the individual decay channels). The errors 
on the numbers of events are statistical (first) and sys- 
tematic (second) where the later is obtained by varying 
the parameterization of the background. 
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FIG. 1: Invariant mass spectra for Knn (n=l,2,3) combina- 
tions in the singly tagged D° decay modes: (a) K + tt~, (b) 
K + tt~tv~tv + and (c) K + ir~ir°. 
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FIG. 2: Invariant mass spectra for mKnir 

(m=0,l,2,n=l,2,3,4) combinations in the singly tagged 
D~ decay modes: (a) K + w~-!r~, (b) K°iy~ , (c) K°R-, (d) 
K+K'TT-, (e) ^ 7r"7r-7r+, (f) K°ty-%°, (g) K + ir--K-ir°, 
(h) K + ir + ir~ ir~ ir~ , (i) ir + ir~ir~ and (j) sum of nine modes. 



C. Doubly Tagged Events 

1. Inclusive decay D -> K*°(K*°)X 

The inclusive decay D — > K*°X is reconstructed on the 
recoil side of the singly tagged D, where the K*° is recon- 
structed through its decay to A~7r + . The regions within 
a ±3cr Md window around the fitted D masses in the in- 
variant mass spectra for mKnir combinations, as shown 
in fig. and fig.|21 are defined as the singly tagged D sig- 
nal regions (D tag region), where <tm d are the standard 
deviations of the mass spectra. Those outside a ±4om d 
window around the fitted D masses are taken as the side- 
band background control regions (D sideband). In the 



estimation of the number of background events in the D 
tag region, the number of the events in the D sideband 
is normalized to the area of the fitted background in the 
D tag region. 

Due to particle misidentification and random combi- 
nation, the invariant mass for A _ 7r + combination could 
enter into mass spectra more than once per event. To 
avoid this problem, only the A"7r + combination with 
the maximum product of the confidence levels for the K 
hypothesis and ir hypothesis is retained in an event. 

Invariant mass spectra for K~ir + combinations ob- 
served on the recoil side of the 5° tags (the mKnn com- 
binations) are shown in fig. for study of the Cabibbo- 
favored decay D° — > K*°X, where (a) shows the mass 
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spectrum for the events with a tagged D°, and (b) shows 
the normalized mass spectrum of the D° sideband events. 
By fitting the invariant mass spectra for K~?t + combi- 
nations with a Gaussian function for the K*° signal and 
a polynomial to describe background, the numbers of the 
observed K*° events are obtained to be 188.5 ± 37.3 and 
92.6 ± 23.5 for the tagged and sideband events, respec- 
tively. In the fit, the mass and width of K*° are fixed at 
0.8961 GeV/c 2 and 0.0507 GeV/c 2 quoted from PDG 0, 
and the detector resolution is set to be 0.0078 GeV/c 2 , 
which is obtained by fitting the invariant mass spectrum 
for K~tt + combinations from both D° and D + decays. 
After subtracting the number of the background events, 
we obtain 95.9 ± 44.1 signal events for the D° — > K*°X 
decay. 

Similarly, fig. 0] shows the invariant mass spectra for 
K~ir + combinations selected on the recoil side of D~ 
tags to study the Cabibbo-favored decay D + — > K*°X. 
Fig. and fig. show the invariant mass spectra for 
K + tt~ combinations on the recoil side of D° and D~ tags 
to study the Cabibbo-suppressed decays D° — > K*°X 
and D + — > K*°X. With the same analysis procedure as 
above, the numbers of the observed K*° /K*° events are 
obtained (Table QJ. After subtracting the number of the 
background events, we obtain 189.1±36.0, 30.8±13.2 and 
12.3 ± 23.3 signal events for D+ -> K*°X, D a -> K*°X 
and D + — > K*°X decays, respectively. 
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FIG. 3: Invariant mass spectra for K tt + combinations ob- 
served on the recoil side of the D° tags for study of the 
Cabibbo-favored decay D° -> K*°X: (a) the mass spectrum 
for the events with a D° tag, (b) the normalized mass spec- 
trum for the D° sideband events. 



2. Efficiencies for D -> K*°(K*°)X 

The efficiencies for reconstruction of the inclusive K*° 
decays of D mesons are estimated with the Monte Carlo 
simulation. The Monte Carlo events are generated as 
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FIG. 4: Invariant mass spectra for K~n combinations ob- 
served on the recoil side of the D~ tags for study of the 
Cabibbo-favored decay D + — > K*°X: (a) the mass spectrum 
for the events with a D~ tag, (b) the normalized mass spec- 
trum for the D~ sideband events. 
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FIG. 5: Invariant mass spectra for K + n~ combinations ob- 
served on the recoil side of the D° tags for study of the 
Cabibbo-suppressed decay D° — > K*°X: (a) the mass spec- 
trum for the events with a D° tag, (b) the normalized mass 
spectrum for the D° sideband events. 



e + e~ — > DD, where D decays into the singly tagged D 
modes and D decays into K*°X. The particle trajecto- 
ries are simulated with the GEANT3 based Monte Carlo 
simulation package of the BESII detector 0- The aver- 
age efficiencies are obtained by weighting the branching 
fractions of D meson decays quoted from PDG and 
the numbers of the singly tagged D events. The efficien- 
cies are 0.1575 ± 0.0015 for the decay D° -* K*°{K*°)X 
and 0.1529 ± 0.0021 for the decay D+ -> K*°(K*°)X. 
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FIG. 6: Invariant mass spectra for K + tv~ combinations ob- 
served on the recoil side of the D~ tags for study of the 
Cabibbo-suppressed decay D + — ► K*°X: (a) the mass spec- 
trum for the events with a D~ tag, (b) the normalized mass 
spectrum for the D~ sideband events. 



TABLE I: Number of K*°/K*° events observed on the re- 
coil side of the D tags, where N and Nb are the number of 
K*° /K*° events observed from the events in which the invari- 
ant masses of the mKn-n combinations are in the D signal 
region and in the D sideband region, respectively, n is the 
number of the signal events for D decays. 



Decay Mode 


N 


N b 


n 


D° -> K*°X 


188.5 ± 37.3 
30.8 ± 13.2 


92.6 ± 23.5 
0.0 ±0.1 


95.9 ±44.1 
30.8 ± 13.2 


D+ -> K* U X 
D+ -> K*°X 


232.5 ± 30.9 
43.7 ± 17.6 


43.4 ± 18.4 
31.4 ± 15.2 


189.1 ±36.0 
12.3 ± 23.3 



IV. RESULTS 

With the numbers of the observed signal events for the 
decay D — ■» K*°X, the numbers of the singly tagged D 
mesons and the reconstruction efficiencies, the branching 
fractions for the Cabibbo-favored decay D — > K*°X are 
determined to be 



and 



BF(D° — > K*°X) = (8.7 ± 4.0 ± 1.2)% (1) 



BF(D+ -> K*°X) = (23.2 ± 4.5 ± 3.0)%. (2) 



These results are consistent with those measured by the 
BES collaboration based on the data taken at ^/s = 4.03 
GeV El. 



For the Cabibbo-suppressed decay D° 
branching fraction is obtained to be 



K*°X the 



BF(D -► K* U X) = (2.8 ± 1.2 ± 0.4)%. 



(3) 



An upper limit on the branching fraction at 90% C.L. for 
the Cabibbo-suppressed decay D + — > K*°X is set to be 



BF(D+ -> K*°X) < 6.6%, 



(4) 



which includes the systematic uncertainty. 

If we treat the observed events as signal events, the 
branching fraction for the Cabibbo-suppressed decay 
D+ -> K*°X is 



BF{D+ -> K*°X) = (1.5l;;g ± 0.2)% 



(5) 



In the measured branching fractions, the first error is 
statistical and second systematic. The systematic error 
arises from the uncertainties in particle identification (~ 
1.0%) 0, in tracking (2.0% per track), in the numbers 
of the singly tagged D° (~ 4.5%) and D~ (~ 3.0%) 0, 
in background parameterization (~ 12%), and in Monte 
Carlo statistics (0.16% for D° and 0.31% for D+). These 
uncertainties are added in quadrature to obtain the total 
systematic error, which is 13.5% for the D° decay and 
13.0% for the D+ decay. 



V. SUMMARY 

Based on an analysis of about 33 pb^ 1 of data col- 
lected with the BES-II detector at the BEPC collider, we 
measured branching fractions for the inclusive Cabibbo- 
favored decay D — > K*°X and Cabibbo-suppressed de- 
cay D -> K*°X. The results_are BF(D° -> K*°X) = 
(8.7±4.0±1.2)%, BF(D+ -> K*°X) = (23.2±4.5±3.0)% 
and BF(D° -> K*°X) = (2.8 ± 1.2 ± 0.4)%. We set an 
upper limit on the branching fraction at 90% C.L. for the 
decay D+ -> K* a X to be BF(D+ -> K*°X) < 6.6%. 
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